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10 Abstract: Forest fires can burn large areas, but can also inject smoke into the upper 

11 troposphere/lower stratosphere (UT/LS), where stakes are even higher for climate, 

12 because emissions tend to have a longer lifetime, and can produce significant regional 

13 and even global climate effects, as is the case with some volcanoes. Large forest fires 

14 are now believed to be more common in summer, especially in the boreal regions, 

15 where pyrocumulus (pvroCu), and occasionally pyrocumulonimbus (pyroCb) clouds 

16 are formed, which can transport emissions into the UT/LS. A major difficulty in de- 

17 veloping realistic fire plume models is the lack of observational data within fire 

18 plumes that resolves structure at a few 100 m scales, which can be used to validate 

19 these models. Here, we report detailed airborne radiation measurements within strong 

20 pyroCu taken over boreal forest fires in Saskatchewan, Canada during the Arctic Re- 
search of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) 
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1 summer field campaign in 2008. We find that the angular distribution of radiance 

2 within the pyroCu is closely related to the diffusion domain in water clouds and can 

' be described by very similar simple cosine functions. We demonstrate with Monte 

4 Carlo simulations that radiation transport in pyroCu is inherently a 3D phenomenon 

5 and must account for particle absorption. However, the simple cosine function prom- 

6 ises to offer an easy solution for climate models. The presence of a prominent smoke 

7 core, defined by strong extinction in the UV, VIS and MIR, suggests that the core 

<8 might be an important pathway for emission transport to the upper troposphere and 

9 lower stratosphere. We speculate that this plume injection core is generated and sus- 

10 tained by complex processes not yet well understood, but not necessarily related di- 

1 1 rectly to the intense fires that originally initiated the plume rise. 
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1 1. Introduction 

2 Forest fires are known to burn vast areas, and also inject smoke into the upper tropo- 

3 sphere/lower stratosphere (UT/'LS), where the stakes are even higher for climate, because 

4 emissions tend to have a longer lifetime, and can cause regional and global climate effects, 

5 as is the case with some volcanic eruptions. Fires of this nature are common in summer in 

6 the boreal regions of Alaska, Canada, and Russia (Fromm et ah, 2010; Soja et ah, 2007; 

7 Stocks et ah, 2003). 

8 Large forest fires release enough energy, both sensible and latent heat, to lift the 

9 smoke-laden air beyond the condensation level, resulting in the formation of cumulus 

1 0 clouds. Due to their association with fires, these clouds are commonly referred to as pyro- 

1 1 cumulus (pyroCu), or pyrocumulonimbus (pyroCb) when convective storms are sufficiently 

12 deep to glaciate and form anvils. The pyroCbs have unique dynamical and microphysical 

13 structures, and constitute an important mechanism for troposphere-stratosphere transport of 

14 biomass burning emissions such as aerosols and trace gases, thereby contributing to raising 

15 their budgets in the UT/'LS (Damoah et ah, 2006; Luderer et ah, 2006; Trentmann et ah, 

16 2006; Fromm et ah, 2005; Fromm and Servranckx, 2003). There is now growing evidence 

17 that wildfires should be considered as additional sources of aerosol and trace gases to the 

18 UT/LS region, which was previously thought to be primarily dominated by volcanic erup- 

1 9 tions (Fromm et ah, 2010; Luderer et ah, 2006). 

20 There is still significant uncertainty on how much of the energy released by combus- 

21 tion contributes to local heating of the atmosphere and is available for convection, and how 

22 much of the energy is lost due to radiative processes. Commonly found estimates for the 

23 radiative energy are between nearly 14% (Wooster et ah, 2005) find 50% (McCarter and 

24 Broido, 1965; Packham, 1969). These estimates are based on laboratory studies or obscrva- 

25 tic ns of small-scale fires, and their application to large-scale crown fires resulting in 
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! pyroCu or pyroCb convection is subject to considerable uncertainty. The interaction of this 

2 radiative energy with atmospheric constituents is also highly uncertain. In the thermal in- 

3 frared, where most of the lire radiation is emitted (Wooster, 2002), aerosols are rather inef- 

4 tlcient absorbers. It is likely that most of the radiative energy from the fire is absorbed by 

5 cloud droplets or gaseous absorption at cloud base or in air masses that are entrained into 

6 the convective plume. This implies that the radiative energy from the fire is trapped in the 

7 lower part of the pyro-convection and therefore contributes to the convective energy. It is 

8 therefore important to consider the radiative processes in detail. 

9 A major challenge in developing realistic fire plume models is the lack of observa- 

10 tional data within forest fire plumes that can be used for validating these models. This is not 

1 1 a big surprise, because making observations within forest fire plumes is challenging at best, 

12 and flying into the “core” of rising smoke can be very risky. The possibility of damaging 

13 expensive instrumentation is also high. Here, we report detailed radiation measurements 

14 from within a pyroCu core taken onboard the NASA P-3B aircraft during the Arctic Re- 

15 search of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) ex- 

16 periment, which took place in Canada in June-July 2008 (Jacob et al., 2010). Such meas- 

17 urements are needed in the investigations of pyro-convection on the spatial scale of indi- 

1 8 vidual events, i.e., tens of meters to a few hundred meters. This would contribute to the 

19 understanding of the vertical and cross-isentropic transport mechanisms responsible for 

20 high-altitude tropospheric and stratospheric injection (Luderer et al., 2006). These radiation 

21 measurements are unique; we are not aware of any previous radiation measurements within 

22 a forest fire plume at the scales described in the next section. 

23 This study focuses on 2 July 2008 and 6 July 2008, the fourth and sixth flights of the 

24 P-3B, respectively, out of Cold Lake, Canada during ARCTAS. These flights provided op- 
portunities for in situ observations of flaming fires, plume evolution and vertical profiles of 
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1 plumes in homogeneous lire outflow regions. Most importantly, they are the only flights 

2 where the P-3B penetrated the pyroCu clouds of active fire plumes. In this study, we em- 

3 phasize data obtained with a multi-wavelength scanning radiometer, NASA’s Cloud Ab- 

4 sorption Radiometer (CAR: Gatebe et al„ 2003: King et al., 1986), which provides unpre- 

5 cedented detail on pyroCu clouds. 

6 The CAR instrument measures scattered light in 14 spectral bands between 0.34 and 

7 2.30 pm and has a unique scanning geometry and is able to scan a scene of the atmosphere 

8 and surface over a 190° scan angle and at high angular resolution (1° instantaneous field of 

9 view). When it is flown inside clouds, it is able to provide a side view in the zenith and na~ 

10 dir directions, and all scattering angles in between from an aircraft platform with as much 

11 as a 5° roll-— assuming no aircraft pitch. Data are sampled simultaneously and continuously 

12 on nine individual detectors. Eight of the data channels for spectral bands from 0.34 to 1.27 

13 pm are always registered during the operation, while the ninth data channel is registered for 

14 signal selected among six spectral channels (1.55-2.30 pm) on a filter wheel. The filter 

15 wheel can either cycle through all six spectral bands at a prescribed interval (usually chan- 

16 ging filter every fifth scan line), or lock onto any one of the six spectral bands, mostly 

1 7 1 .656, 2. 103, or 2.205 pm, and sample it continuously 

18 2. Observational Assessment of the PyroCu 

19 Some of the biggest Saskatchewan fires during the 2008 summer ARCTAS campaign 

20 burned west of Mcintosh Lake and east of Lake Athabaska in Canada. On 2 July (P-3 Flight 

21 # 1 8/ CAR Flight # 20 1 7) and 6 July ( P-3 Flight #20/CAR Flight # 20 1 9) the conditions for 

22 fires seemed right, as a cold front moved over the region setting up conditions conducive to 

23 sending hot smoke plumes billowing into the upper atmosphere. The P-3B aircraft, with its 

24 suite of instruments, was in the air and managed to acquire very good fire data. Below, we 
provide a general description of the 2 July flight and the primary instruments that were op- 
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! crating during the flight. We will not describe the 6 July flight because of similarities be- 

2 tween the two flights, however, full details for this flight can be found on the CAR website 

3 (http://car.gsfc.nasa.gov/data/index.php7id~! 1 5&mis_id=8&n=ARCTAS). Through these 

4 descriptions, we attempt to paint a picture of the nature and character of the fires and in ad- 

5 dition, create a preamble for further work on pyroCu and their effects on climate. 

6 2.1 General Description of Flight 201 7 , 2 July 2008 

7 Our objectives for this flight included flying a segment under the Terra satellite, co- 

8 ordinated with the NASA B-200 carrying the High Spectral Resolution Lidar (HSRL) and 

9 the Research Scanning Polarimeter (RSP) instrument) with radiation measurements under 

10 clear sky conditions stacked at different altitudes, and in situ observations of flaming fires, 

1 l showing plume evolution and vertical profiles in homogeneous fire outflow regions. The 

12 flight took 6 hours and 20 minutes, with take off at 15:3 1 UTC (9:31 a.m. local time; see 

13 flight log: http:/7car.gsfc.nasa.gov/data / index.php?id=113&mis_id= : 8&n=ARCTAS). Fig. 

14 la shows the flight path from Cold Lake, heading northeast towards Reindeer Lake. We 

15 made measurements collocated with the B-200 for about 50 minutes beginning at 16:25 

16 UTC, and observed active fires of different intensities under scattered cloud conditions. We 

17 penetrated an active fire plume near Mcintosh Lake at two altitudes, 1303 m (19:18 UTC) 

1 8 and 294 m (19:49 UTC) above the local surface, and flew along the plume axis for about 80 

19 km to study the evolution of aerosol microphysics during downwind transport. Between the 

20 two penetrations, we obtained radiation measurements in a vertical profile in the most 

21 homogeneous part of the plume and flew across-plume at a low-level. We also sampled 

22 outflow from a pyroCu cloud and penetrated the cloud to sample the inside. Fig. lb shows 

23 an image taken by the Moderate Resolution Imaging Spectroradiometer (MODIS) flying 

24 aboard the Aqua satellite at 19:35 UTC over the region highlighted in the white box in Fig. 

25 la. The cloud fields seen in the satellite image are most likely associated with the active 
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1 fires in the region. Note that the green squares represent fire locations detected by MODIS 

2 on the Terra satellite, taken during the overpass earlier that day, at 17:55 UTC. The red 

3 squares represent fire locations based on high radiance, relative to background, as seen with 

4 the MODIS/ Aqua 4-pm channel during the overpass (e.g. Giglio, 2003). 

5 Figs. 1c and Id show CAR quicklook images selected from a section of the flight 

6 where the aircraft went through a pyroCu core during the flight of 2 July (the green line 

7 highlights the flight path through the pyroCu, whereas the orange dot in Fig. lb marks the 

8 location of the plume core). A combination of different CAR bands representing red, 

9 green, and blue colors helps to differentiate between clouds and smoke and to see through 

10 smoke, revealing land surface features such as green vegetation and dark water bodies (see 

11 Fig. lc: false-color RGB(0.47, 0.38, 0.34) pm; Fig Id: false-color RGB( 1 .66, 1.04, 0.38) 

12 pm). As shown in Fig. Id, the fire was widespread, covering hundreds of hectares, and it 

13 took the aircraft - 15 seconds to cross the core (the core covers a ground area ~2 km, where 

14 the measured aircraft groundspeed varied between 127 ms -1 and 138 ms -1 ). The aircraft 

15 altitude changed dramatically while traversing the core, rising by 93 m, from 2747 m to 

16 2840 m above mean sea level, an indication of strong updrafts >8 ms' 1 . 

17 The height of the 2 July PyroCu top above ground was > 3000 m, and the boundary 

1 8 layer (BL) height at the fire location and time of aircraft data acquisition was 1976 m. The 

19 BL was obtained by interpolating the values for 12 noon and 3 pm local times from the 

20 NASA Goddard Earth Observing System (GEOS-5) data provided by the Global Modelling 

21 and Assimilation Office (GMAO). Therefore, the pyroCu extended well above the boun- 

22 dary layer. Recent studies based on digitized plume height analysis from the Multi-angle 

23 imaging SpectroRadiometer (M1SR) instrument aboard the Terra satellite have shown that 

24 of order 10% of the plumes are injected higher than 0,5 km above the BL; out of eight 

25 North American biome types evaluated, fire plumes from boreal forests (which include the 
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1 study area) are typically injected the highest, with a median injection top height of -1200 m 

2 (as compared to -800 m for croplands; Kahn et al., 2008; Val Martin et al., 2010). Those 

3 studies also found a correlation between the plume injection height and lire radiative power 

4 (FRP) measurements from MODIS (Ichoku et al, 2008), indicating that in addition to at- 

5 mospheric stability structure, tire intensity (represented by FRP) has a strong influence on 

6 plume injection height; plume top heights that reached the free troposphere (FT) were pro- 

7 duced mainly by high intensity fires having a median FRP of order 500 MW. Table 1 

8 shows the FRP values of the MODIS fire pixels (depicting their relative intensities) from 

9 Aqua coinciding with 2 July 2008 flights (in space but not in time). The FRP data record 

10 corresponding spatially to the pyroCu plume penetrated by P-3B is shown in boldface char- 

1 1 actors. However, because of the large difference in time between the FRP measurement 

12 from Aqua-MODIS (19:35 UTC) and P-3B penetration of the pyroCu (-21:05 UTC), this 

1 3 FRP value is much less than what would have been required to generate such a high plume 

14 and the associated pyroCu. It is probable that the fire became much stronger by the time of 

15 the pyroCu measurement. For the P-3B pyroCu penetration on 6 July 2008, no fire pixels 

16 were detected by MODIS within a few km of the pyroCu. It is believed that the pyroCu 

17 may have moved substantially away from the fire at the lime of its penetration, making it 

1 8 impossible to link it to any specific fire detected by MODIS. Therefore, to establish a cause 

19 and effect relationship of such transient features as fires and their associated plumes and 

20 pyroCu, it is crucial to conduct airborne measurements in synchronicity with satellite over- 

2 1 passes and ground-based measurements whenever possible. 

22 2. 2 NASA P-3B Payload during A RCTAS 2008 

23 Fig. 2a shows the platform we used in this campaign, the NASA P-3B aircraft. The 

24 picture was taken at NASA Ames on 24 June 2008 just before a transit flight to Cold Lake, 
Canada at the start of the A RCTAS summer campaign. The P-3B payload included 10 pri- 
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! snary instruments (see Table 2) for measuring aerosol optical depth, aerosol extinction and 

2 scattering, aerosol size distribution (Russell et al., 1999; Clarke et al., 2007), spectral (0.380 

3 - 2.200 pm) and broadband downwelling and upweiling solar (0.2-3. 6 pm), and 1R irradi- 

4 ance (4.5-42 pm; Pilewskie et al., 2003; Rucholtz personal communication: 

5 http://airborncscience.nasa.gov/instrument/BBR), angular distribution of scattered radiation 

6 in different directions (Gatebe et al., 2003), cloud condensation nuclei (Roberts and Nenes, 

7 2005; Lance et ah, 2006), and NCR and carbon monoxide (Stravva et al., 2007; Provencal et 

8 al. 2005). The base meteorological and navigational measurements were provided by the 

9 Project Data System, which is the primary navigational data system for P-3B investigators. 

10 Fig. 2b shows a schematic of the CAR instrument (Gatebe et al., 2003; King et al., 

1 1 1986), highlighting its main features. Fig. 2c shows a pyroCu taken onboard the NASA P- 

12 3B on 2 July 2008 near Reindeer Lake, Saskatchewan, Canada. The aircraft Hew around 

13 these fires, sampling dean air upwind and the polluted air downwind. It then flew through 

14 the smoke at various altitudes, probing and charting changes in concentration and the 

15 "chemical evolution" of the compounds emitted by the fires. The plane also flew into the 

16 i4 core” of the rising fire plumes, where the CAR instrument managed to acquire good data 

17 as discussed in this study. These represent the first such measurements in the core of a 

18 pyroCb of which we are aware. (See pictures from the 2 July flight on the CAR website at 

19 http://car.gsfc.nasa.gov/data/photos.php7mis id=8&n=ARCTAS&d=2008-07- 

20 02&f=20 1 7&l=h and from the 6 July flight at 

21 http://car.gsfc.nasa.gov/data/photos.php?mis_id=8&n=ARCTAS&d=2008-07- 

22 06&f--2019&l=h. 

23 In the next section, we will present an analysis of the CAR measurements inside 

24 pyroCu and results from a Monte Carlo simulation of the radiation field within optically 

dense media analogous to those observed above. 
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1 3. Results and Discussion 

2 3.1 Analysis of CAR Measurements 

3 The core of the pyroCu as seen in Figs, lc and Id stands out as a very dark region of 

4 the plume, where radiation seems to be uniformly distributed and dominated by strong ab- 

5 sorption in all directions. The CAR instrument measured very low spectral radiance values 

6 (< 1 Wm'^m^sr 1 in the ultraviolet and visible bands having X < 0.5 pm, and at X = 1.6 

7 pm, and 2-5 Wm" 2 pm‘*sr l for bands between 0.5 pm and 1.3 pm), independent of direc- 

8 tion. The radiance distributions within the core at two zenith angles: 5° and 175° at 0.472 

9 pm, shown in Fig. 3 demonstrate that there is hardly any difference between the downwel- 

10 ling and upwelling radiance within the core region. The a, b and c images in this figures 

1 1 were taken by a forward camera aboard the NASA P-3B at the entry, core, and exit points 

12 in the pyroCu. Note the brownish color of the plume inside the core region. The radiance 

13 decreases by about two orders of magnitude from outside to the central part of the core, ir- 

14 respective of the direction and wavelength. We speculate that the core region plays an im- 

1 5 portant role as the main physical pathway through which emissions from fires are trans- 

16 ported into the upper troposphere/lower stratosphere (UT/LS), but this observation needs to 

17 be validated. 

18 Fig. 4 compares relative spectral radiance as a function of zenith angle for internally 

19 scattered radiation measured deep inside an optically thick water cloud, pyroCu and a thick 

20 “pure” smoke, which is defined by lack of reflectance sensitivity in the near-infrared bands 

21 (X > 0.87 pm), at selected CAR wavelengths. According to King et ah (1990), radiance 

22 measurements from an individual scan of the CAR are considered to be in the diffusion 

23 domain, if three conditions are satisfied: (i) the zenith radiance exceeds the nadir radiance, 

24 (ii) the maximum deviation from the theoretical cosine curve is less than or equal to 5% of 
the mean amplitude, and {iii) the number of times the deviations from a cosine distribution 
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! changes sign is greater than or equal to 4, which helps to assure that fluctuations are ran- 

2 dom, and not systematic drifts. These criteria pertain to the diffusion domain of an optically 

3 thick media if the scattering is conservative (no absorption), and thus are to be expected for 

4 liquid water clouds at wavelengths around 0.67 pm (as seen in Fig. 4a). The theoretical 

5 curve is generated from measured zenith and nadir radiances using the following cosine 

6 relationship 

7 7(r,cos0) =a + bco$0 (1) 

8 where 

9 a = ~(/ 0 (t,1) + / 180 (t,-1)), 

10 b = — (/ 0 (r,l) -/ 18 o(t,-1)) (la) 

1 1 and / 0 and I m are the zenith and nadir radiances at 0 = 0° and 0 = 1 80°, respectively. From 

12 Eq. la, we came up with a general equation for any two downward and upward viewing 

1 3 angles in the form 

14 a = ^(/„.(t,cos0') + cos 0 ')), 

15 b = —^—(I e iT,cosd')-r im _ ( ,iT,-cosd')) (lb) 

2 cos O’ 

16 As such, Eq. lb is not restricted to using zenith (0 = 8’ = 0°) and nadir (0 = O’ = 180°) in- 

17 put radiances, unlike Eq. la, and can be used to describe the radiation Held in optically 

1 8 thick media given any two downward and upward viewing directions, or 8’ = 90° (where b 

19 = Ic)q(t) ) and any other direction. This general formula can be applied to CAR observa- 

20 tions away from the zenith and nadir directions, and therefore can accommodate the reali- 

2 1 ties of data acquisition on an airborne platform that has a nonzero pitch and/or roll angle. 

22 For the liquid water cloud case, obtained on 29 June 2008 (20:00 UTC) in Canada, 

23 the angular radiance field at the shortest wavelengths can be described by a cosine function 

24 expected for conservative scattering in the diffusion domain (King, 1981; King et at., 

25 1990), where the cloud is expected to have a small to negligible amount of absorption. This 
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! is demonstrated in Fig. 5 (liquid water cloud case), where the angular distribution of the 

2 radiation field obtained from measurements (continuous curve, abstracted from Fig. 4a) is 

3 compared to that expected in the diffusion domain for conservative scattering (dotted 

4 curve) at X = 0.472 pm. We note that the angular radiance field becomes increasingly ani- 

5 sotropic as absorption increases in the CAR bands at wavelengths exceeding 1 pm (Fig. 

6 4a), because stronger absorption implies that a smaller portion of do vvn welling solar radi- 

7 ation can be reilccted back up from below. In this case, the three conditions for the diffu- 

8 sion domain at a conservative scattering wavelength were met. King et ah, (1990) show that 

9 when the three conditions are assured, quantitative information about cloud absorption 

10 properties from the angular distribution of scattered radiation can be derived, dependent 

1 1 only on the optical properties of the medium and possibly on the reflectivity of the under- 

1 2 lying surface, but independent of the solar zenith angle and solar irradiance. 

1 3 The problem now is to describe the radiation pattern inside an optically thick pyroCu, 

14 the core of a pyrcoCu, and smoke, which seem to be increasingly anisotropic as absorption 

1 5 tends to dominate at the shorter wavelength and is not negligible. However, the diffusion 

16 domain theory for optically thick clouds with modified a and b (Eq. lb; Fig 4b), can also 

17 describe the angular radiance distribution in a pyroCu as shown in Fig. 5 (blue solid curve 

18 represents measurements, green dotted line represents the fit), but b must assume a negative 

1 9 value, because in this case, the relative radiance is highest at nadir rather than at zenith. The 

20 cosine function can also represent the pyroCu core case (Fig. 4c), but an additional sin0 

21 term is needed (Eq. 2) to reproduce the bow shaped curves in Figs. 4c and 5 (note that the 

22 solid lines represent observations, while dotted lines represent fit to the observations). 

/(r,cos 0) = (a + i?cos0) + (c + dsin0) (2) 

which requires only the knowledge of the zenith angle, 0, and some constants. The param- 
eters a and b are easily determined from Eq. la & lb. The last two terms represent a math- 
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! ematical ansatz, where the term c is obtained from zenith and nadir radiances following Eq. 

2 lb, and d was determined through trial and error until the best fit was found (d = -0.1). It 

3 is interesting to see that radiation transport in an optically thick pyroCu core can be de~ 

4 scribed by simple geometrical functions, whereby the first two terms represent scattering 

5 processes, while the third and fourth terms represent absorption processes. To gain better 

6 understanding of radiation transport in a pyroCu, we next attempt to reproduce the radiation 

7 pattern observed by the CAR instrument using ID and 3D Monte Carlo simulations. 

8 3.2 Monte Carlo Simulation 

9 We now explore the optical characteristics of pyroCu using I D and 3D Monte Carlo 

10 radiative transfer simulations. The Monte Carlo model is an appropriately modified version 

1 1 of the code used in several earlier studies, including Vamai and Marshak (2001), and was 

12 also tested in the Intercomparison of 3D Radiation Codes (I3RC) project (Cahalan et al.. 

13 2005). The actual code used in this paper was also tested through comparisons with 

14 DISORT (Discrete Ordinates Radiative Transfer) model (Stamnes et al. 1988) runs for ab- 

15 sorbing, plane-parallel clouds. To keep simulation uncertainties below 0.5%, each result is 

16 based on 4 million photon trajectories in backward Monte Carlo simulations. 

1 7 Because the purpose of these Monte Carlo simulations is to understand the basic be- 

18 haviours in Fig. 4 (and not to simulate the exact values observed by CAR), the simulations 

19 use a simple setup (Fig. 6). The calculations assume a 2D slab cloud that is infinite in the 

20 cross-sun (Y) direction, lies between 1 km and 2 km altitudes, and is l km wide (AZ = AX 

21 =1 km). The extinction coefficient is 100 km" 1 , so the cloud optical thickness is 100 in 

22 both Z and X directions. Because we simulate radiances deep inside the cloud, where 

23 multiple scattering minimizes the details of the scattering phase function, we assume the 

24 cloud particles have a Flenyey-Greenstein phase function with an asymmetry parameter of g 

25 = 0.85. The solar zenith angle is taken as 30 ' and the surface albedo as 0.2. For simplicity. 
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1 no atmospheric effects — Rayleigh scattering or gaseous absorption — are considered. 

2 Fig. 7 shows simulation results that can shed light onto the behaviours in Figs. 4a. 4b. 

3 and 4c. (At this stage we do not attempt to simulate the complex behaviours in Fig. 4d.) 

4 The black dashed line displays a qualitatively similar behaviour to the clean cloud in Fig. 

5 4a and shows the general behaviour of ID radiation calculations, which assume horizon- 

6 tally homogeneous and infinite (plane-parallel) clouds. The line illustrates that in ID cases, 

7 the radiance of the highly diffuse radiation deep inside clouds is always greater looking up 

8 than looking down regardless of the amount of particle absorption. This is because the con- 

9 servation of energy and the lack of light sources below imply that the upwelling Dux can 

10 never exceed the downwelling flux. This particular simulation was based on a single scat- 

1 t tering albedo of 0.975, representing stronger absorption than typically occur in liquid water 

12 clouds in the visible near-infrared (cf. Fig. 4a), but with a shape typical of what occurs in 

1 3 absorbing liquid water clouds. 

14 The other three curves in Fig. 7 are from 3D simulations. Energy conservation re- 

15 quires the scene average net radiation to How downward even in 3D situations. However, 

16 this trend can reverse locally inside clouds in 3D situations because large amounts of radi- 

17 ation can reach the surface easily in clear areas between clouds and, after reflection from 

18 the surface, can enter through the cloud base and create an upwelling flow of radiation in- 

19 side the cloud. Jn addition, sunlight entering through a sunlit cloud side can create a hori- 

20 zontal net flow toward the shadowy side (Fig. 6). 

2 1 The blue squares in Fig. 7 show that the downward flow entering through the cloud 

22 top still dominates our 3D case if there is no absorption. The red circles, however, show- 

23 that absorption can drastically change the radiation pattern by depleting much of the 

24 downwelling radiation wending its way through the cloud. As a result, the upwelling flow- 

25 of radiation reflected from the surface can match the weakened downwelling flow, creating 
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1 comparable intensity peaks looking up (0= 0°) and down (0=180°). This is similar to the 

2 behaviour in Fig. 4c (pyroCu core). Lastly, the green dashed curve in Fig. 7 shows that for 

3 a lower altitude cloud absorption reduces the downweliing (upwelling) radiation more 

4 (less), and so upwelling radiation can fully dominate and create a trend similar to that in 

5 Figure 4b (pyroCu). 

6 Finally, Fig. 8 more explicitly explores the role of absorption and scattering proper- 

7 ties in shaping radiation patterns. Fig. 8a shows that even relatively weak absorption (SSA < 

8 0.99) can allow upwelling radiation to match or even exceed the downweliing flow in a 3D 

9 situation. In turn. Fig. 8b shows that increasing the number of particles (i.e., the extinction 

10 coefficient) yields similar effects to increasing their absorptivity, as both changes make it 

1 1 harder for radiation to reach an instrument deep inside the cloud. This implies that radiation 

12 patterns similar to those observed in pyroCu and pyroCu core (Figs. 4b and 4c) can arise in 

13 thick clouds even if absorption is relatively weak, although stronger absorption is required 

14 in thin clouds. The interplay between scattering and absorption can also allow upwelling 

15 radiation to dominate even in the upper parts of clouds if particles are distributed unevenly. 

16 For example, if the cloud extinction coefficient increases with altitude (e.g., Raga and 

17 Jonas, 1993) — e.g., as ascending droplets grow through condensation— strong absorption or 

18 scattering above can greatly reduce the downweliing radiation and allow the upwelling ra- 

19 diation reflected from the ground to dominate even at high altitudes. 

20 4. Conclusions 

2 1 This study provides an analysis of airborne radiation measurements with the CAR in- 

22 strum ent within strong pyroCu over boreal forest fires in Saskatchewan, Canada during the 

23 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARC- 
TAS) summer field campaign in 2008. 


24 
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1 We found the pyroCu core to be very interesting having very low light levels from the 

2 UV to near IR regions of the electromagnetic spectrum, and very uniform radiance field in 

3 all directions. Spectral radiances measured within the core were very small < 1 Wm‘ 2 pm" 

4 •sr 1 in the ultraviolet and visible bands, X < 0.5 pm, and at a,= 1.6 pm, but substantially lar- 

5 ger (2-5 Wnr 2 pnr'sr l ) for 0.5 pm < X < 1.3 pm, independent of direction. We found the 

6 angular distribution of radiance within the pyroCu to be related to the diffusion domain in 

7 water clouds that can be described by very similar cosine functions with an additional sine 

8 term. We demonstrate with Monte Carlo simulation that radiation transport in pyroCu is 

9 inherently a 3D problem and must account for strong particle absorption. The proposed 

10 simple description of the radiation distribution within the pyroCu or pyroCu core may lead 

1 1 to better understanding of these cloud systems and provide a simpler solution for their in- 

12 elusion in climate models. Based on the low radiation values measured in the core of the 

13 pyroCu, we speculate that the formation of the dark core is an important feature for estab- 

1 4 lishment of a possible pathway for transportation of fire emissions to the upper troposphere 

15 and lower stratosphere, but this remains to be validated. 
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1 Table 1 

2 Locations and FRP values of lire pixels detected by Aqua-MODIS on 2 July 2008 at 19:35 

3 UTC in the vicinity of the measured pyroCu, as indicated in Fig. lb 


Latitude 

(°N) 

Longitude 

(°W) 

Fire Radiative Power 
(MW) 

55.685 

105.283 

108.273 

55.682 

105.299 

112.081 

55.677 

105.331 

37.455 

55.688 

105.320 

84.818 

55.706 

105.320 

105.710 

55.733 

105.334 

20.674 

55.778 

105.235 

24.367 

55.773 

105.267 

18.737 

55.803 

105.197 

84.151 

55.838 

105.150 

25.296 

55.835 

105.166 

17.619 
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Table 2 

P-3B Aircraft Instrumentation for ARC T' AS 


Instrument Name Dataproducts 

Technique 

Principal 

Investigator 


AER03X Aerosol optical prop- 

Cavity ring-down ex- 

A. Strawa. NASA 

Extinction and scattering - 2 a 

erties, N0 2 mixing 

tinction; Sample rate: 

ARC 

sensitivity ~ 0.2 Mnr 1 . Extinction 

ratio 

2 sec. (Strawa et al, 


and scattering - 2 a precision 


2007). 


0.1 VI nr' 1 . N0 2 mixing ratio - 2 a 


sensitivity better than ! ppb 


Ames Airborne 

Aerosol optical depth 

Tracking Sun photom- 

J. Redemann. 

Slant OD -0.002 (±0.01) 

Tracking Sunpho- 

and extinction. 

eter. 0.354-2. 138 pm 

BAERI/NASA 

Slant WV -0.0005 to 0.006 gem' 2 

n '-meter (A ATS) 

water vapor column 
and profile 

(Russell et al. 1990) 

ARC 

(±8%) 

Broad Band Radi- 

Radiation fluxes (so- 

Radiometer 

A. Bucholtz, 

Estimated Accuracy: 3-5% 

ometers (BBR) 

lar: 0.2-3. 6 pm & 
IR: 4.5-42 pm) 


NRL 


Cloud Absorption 

Radiance, aerosol. 

Angular (cross-track 

C. Gatebe. 

Radiance absolute accuracy <5% 

Radiometer 

cloud & Earth sur- 

zenith to nadir) scan- 

USRA/NASA 


(CAR) 

face properties 

ning radiometer, 
0.340-2.301 pm 

GSFC 


Continuous-Flow 

Cloud Condensation 

Supersaturation gener- 

A. Nenes. Georgia 

SS Range: 0.20 - 0.60%. ±0.05% 

Streamwise Ther- 

Nuclei (CCN) 

ated by relative diffu- 

Tech 


mal Gradient 


sion of water vapour 



Cloud Condensa- 


and heat (Roberts and 



tion Nuclei 


Nenes. 2005) 
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(CCN) Counter 
(CFSTGC) 
Carbon monoxide 

CO mixing ratio 

off-axis integrated 

J. Podolske, 

By Attenuated 

(mole fraction) at a 

cavity output spec- 

NASA/ARC 

Laser Transmis- 

I -Hz rate based on 

troscopy - TLAS 


sion (COBALT) 

measured absorp- 

(Provencal et al. 


Hawaii Group for 

tion, gas tempera- 
ture. and pressure 
using Beer's Law 
aerosols (Aerosol 

2005) 

OPC. PSAP. TDMA. 

A. Clarke. U. 

Environmental 

number, size, com- 

CN counter. ToF-AMS. 

Hawaii 

Aerosol Research 

position. volatility 7 . 

SP2. nephelometer 


(HiGEAR) 

Project Data Sys- 

optical properties; 
Clarke et a!.. 2007; 
Shinozuka et ah. 
2011, Howell et ah, 
2006) 

base meteorological 

See reference: Barrick 

John Barrick 

tem (PDS) 

(P. T. RH) and 


NASA/LaRC 

Solar Spectral Flux 

navigational meas- 
urements 
Solar spectral flux 

Spectrometer (380- 

S. Schmidt. U. 

Radiometer 


1 700 nm) with nadir 

Colorado 

(SSFR) 


and zenith hemispheric 
collectors 



Precision: 0.2 ppbv (1 -s averaging 
time), uncertainty <1 .0% 


See references: Shinozuka et ah, 
( 2011 ). 


See references: Barrick 


Absolute accuracy 3-5%. 
Precision 1% 
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Fig. 1. (a) NASA P-3B flight track 2 July 2008 during the ARCTAS 2008 summer 

campaign, (b) MODIS/Aqua image taken on 2 July 2008 over the area 
marked by the rectangle in (a). The red (green) squares are fire locations 
during Terra (Aqua) overpass at 19:35 UTC (17:55 UTC), and the orange 
square marks the location of the fire plume that was penetrated by the P- 
3B, (c)-(d) CAR quick-look images of the fires taken aboard the NASA P- 
3B over the flight section highlighted in green in (b). A combination of 
red, green, and blue colors using different CAR bands helps differentiate 
cloud from smoke, and to see through smoke and identify land surface 
features such as green vegetation and dark water bodies. The red, green, 
and blue color assignments in (c) are 0.47, 0.38, 0.34 qm and for (d) 1.66, 
1.04, 0.38 uni. 



(a) NASA P-38 Aircraft 


(c) Pyrocumuius Clouds Studied 
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Fig. 2. (a). The NASA P-3B at NASA Ames Research Center, California, USA in 

June 2008 during the ARCTAS field experiment, (b) Schematic of NASA's 
Cloud Absorption Radiometer (CAR), which is mounted in the nose of 
the NASA P-3B aircraft, (c) Picture of the fire plume that the P-3 pene- 
trated on 2 July during ARCTAS. (d) Specifications for the CAR, which 
contains 14 narrow spectral bands between 0.34 and 2.30 pm. The CAR 
flew nine fire missions aboard P-3B in Canada between 28 June and 10 
July 2008. 


(b) CAR Schematic 


(d) Cloud Absorption Radiometer (CARS Parameters 


Angular Scan Range 
Instantaneous field of view 
Pixels per scan line 
Scan rate 

Spectral channels (pm) 
bandwidth (FWHM) 


190“ 

17.5 mrad (1°) 

382 

1 .67 lines per second ( 1 00 rpm) 

14 * (8 continuously sampled and last 
six In filter wheel) 

0.340(0.009), 0.381(0.006), 0.473(0.021), 
0.683(0.021), 0.871(0.022), 1.037(0.021), 
1.222(0.023), 1.275(0.023), 1.564(0.032), 
1.657(0.042), 1.738(0.040), 2.105(0.045), 
2.202(0.043), 2.303(0.044) 




Time (UTC) 


from a forward camera aboard the NASA P-3B, and coincide with entry, 
core penetration, and exit points of the pyroCu, respectively. The core 
was very dark, and brownish in color. The distance from the core to the 
edge of the pyroCu is about 318 m (assuming average aircraft speed of 
132.5 m/s). 




Fig. 4. Relative spectral intensities from CAR as a function of zenith angle for 
internally scattered radiation measured inside (a) liquid water cloud, (b) 
pyroCu, (c) core of pyroCu, and (d) thick smoke from smoldering fires 
with no cloud present. The UV channels are not included in (b) and (c) 

because of low values of signal-to-noise values. The effects of surface 
reflectance are clearly seen in the "pure smoke" case (d) at wavelengths, X, 
>0.68 ,um. 















0 50 100 150 200 

Scan Angle (°) 

Fig. 5. Observed and fitted relative intensities at 0.472 jUm represented by solid 
and dashed lines, respectively, for water cloud, pyroCu and the pyroCu 
core. The fitted lines for water cloud and pyroCu were generated by a 
cosine function given in Eq. lb, with parameter "b" assuming a negative 
value for the pyroCu case. The pyroCu core was fitted by Eq. 2. 




Illustration of solar illumination paths responsible for the observed radia- 
tion patterns in Figs. 4a-c. 


Fig. 6. 
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Fig. 7. Simulated CAR data used for understanding the observed be- 
haviors in Fig. 4. The figure shows simulated CAR data for an azimuth 
perpendicular to the sun (90° relative azimuth). The black, red, and blue 
curves are for the horizontal center of the cloud at 1.4 km altitude. The 
green curve, simulating behaviors observed at a different fpyroCu) seg- 
ment of the flight path, is for a point 200 m away from the sunlit cloud 
edge at 1.3 km altitude. 
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Fig. 8. Impact of (a) variations in single scattering albedo (SSA) at a constant 
cloud extinction coefficient (100 km -1 ); and (b) variations in cloud extinc- 
tion coefficient at a constant SSA (0.97). The simulations are for 85° rela- 
tive azimuth at a point 300 m away from the sunlit cloud edge at 1.3 km 
altitude. 




